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ABSTRACT 

We derive, adopting a direct method, the luminosity function and the formation rate of long Gamma Ray Bursts through a complete, 
flux-limited, sample of Swift bursts which has a high level of completeness in redshift z (~ 82%). We parametrise the redshift evolution 
of the GRB luminosity as L = Lo(l +z)* and we derive k = 2.5, consistently with recent estimates. The de-evolved luminosity function 
<p(Lo) of GRBs can be represented by a broken power law with slopes a = -1.32 ± 0.21 and b = -1.84 ± 0.24 below and above, 
respectively, a characteristic break luminosity Lo.b = 10^* ''^*°'*^ erg/s. Under the hypothesis of luminosity evolution we find that the 
GRB formation rate increases with redshift up to z ~ 2, where it peaks, and then decreases in agreement with the shape of the cosmic 
star formation rate. We test the direct method through numerical simulations and we show that if it is applied to incomplete (both in 
redshift and/or flux) GRB samples it can misleadingly result in an excess of the GRB formation rate at low redshifts. 
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1. Introduction 

Since the discovery of Gamma Ray Bursts (GRBs), one of the 
most important questions was related to their distance scale (i.e. 
whether galactic or cosmological) which had immediate impli¬ 
cations on their associated luminosities and energetics. Through 
the afterglow detection (Costa et al. 1997; Van Paradijs 1997) 
and first redshift measurements, GRBs were proven to be cos¬ 
mological sources with large isotropic equivalent luminosities 
exceeding, in few cases, 10^"^ erg s '. The pinpointing of the 
GRB afterglow, made available by the fast slewing of the Swift 
satellite (Gehrels et al. 2004), coupled with intense efforts to ac¬ 
quire early time optical spectra from ground, allowed us to mea¬ 
sure the redshifts z of GRBs with an average efficiency of 30%. 
Among these, GRB 090423 (with a spectroscopic z — 8.2 - Sal- 
vaterra et al. 2009a; Tanvir et al. 2009) and GRB 090429B (with 
photometric redshift z = 9.4 - Cucchiara et al. 2011) represent 
the farthest objects of stellar origin known so far. 

Two of the key properties characterising the population 
of GRBs are (a) their cosmic rate fiz) (GRB formation rate, 
GRBFR hereafter), representing the number of bursts per unit 
comoving volume and time as a function of redshift, and (b) 
their luminosity function 0(L) (LF hereafter), representing the 
relative fraction of bursts with a certain luminosity. Here, with 
(p{L) we refer to the differential luminosity function defined as 
dN{L)ldL. 

Recovering i/^(z) and (p{L) of GRBs allows us to test the na¬ 
ture of their progenitor (e.g. through the comparison with the 
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cosmic star formation rate) and to study the possible presence of 
sub-classes of GRBs at the low end of the luminosity function 
(e.g. Liang et al. 2007, see also Pescalli et al. 2015). These two 
functions have been derived for the population of long GRBs 
(e.g. Daigne et al. 2006; Guetta & Della Valle 2007; Firmani et 
al. 2004; Salvaterra & Chincarini 2007; Salvaterra et al. 2009b, 
2012; Wanderman & Piran 2010; Yu et al. 2015; Petrosian et al. 
2015) through different methods and samples of bursts (§2). For 
the population of short GRBs, if{z) and (p(L) have been less se¬ 
curely constrained (e.g. Guetta & Piran 2005,2006; Nakar2006; 
Berger et al. 2014; D’Avanzo et al. 2015) due to the limited num¬ 
ber of bursts with measured redshifts. 


However, fiz) and f(L) cannot be derived straightforwardly 
using all GRBs with known redshift since these samples are af¬ 
fected by observational biases. Specific methods that correct for 
such biases should be adopted. The main approaches that have 
been used so far (§2) agree on the shape of the luminosity func¬ 
tion (typically represented by a broken power law) but lead to re¬ 
markably different results on the cosmic GRB rate (particularly 
at low redshifts). Independently from the method used to recover 
these two functions, most of the previous studies (see however 
Salvaterra et al. 2012) adopted either heterogeneous samples (i.e. 
including GRBs detected by different satellites/GRB detectors 
which have different sensitivities) and/or incomplete samples. 
In particular, incompleteness is induced by several effects such 
as the variation of the trigger efficiency and the redshift mea¬ 
surement. Accounting for such instrumental effects is extremely 
difficult in practice. 
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An alternative is to work with complete samples, at the ex¬ 
pense of the number of GRBs in the sample. Salvaterra et al. 
(2012) (SI2) defined a complete flux-limited sample of GRBs 
(called BAT6) detected by Swift which, despite containing a rel¬ 
atively small number of GRBs, has a high redshift completeness 
and has been extensively used to test various prompt and after¬ 
glow properties of GRBs in an unbiased way (Campana et al. 
2012, D’Avanzo et al. 2012, Ghirlanda et al. 2012, Melandri et 
al. 2012, Nava et al. 2012, Covino et al. 2013, Melandri et al. 
2014, Vergani et al. 2014). 

The aim of this work is to derive 0(L) and f{z) through this 
complete sample of GRBs detected by Swift. We summarise the 
main different approaches that have been used in the literature to 
derive the luminosity function and the formation rate of GRBs 
(§2) and present the updated sample used in this work in §3. We 
adopt the C direct method (Lynden-Bell et al. 1971) to derive 
the (p(L) and if{z) and compare it with previous results in §4. 
Throughout the paper we assume a standard ACDM cosmolog¬ 
ical model with Q,„ = 0.3 and Qa = 0.7 with Hq - 70 km s ' 
Mpc '. We use the symbol L to indicate the isotropic equivalent 
luminosity Ljso omitting for simplicity the subscript “iso". 


2. 0(L) and of long GRBs 

The number of GRBs detectable by a given instrument above its 
sensitivity flux limit S can be expressed as: 


m>s)^ 


Q.T 

An 


S) 

«-^0 L \ 


f(L, z) - -- — dL dz 

(s,z) (I + z) dz 


( 1 ) 


where G and T are the instrument held of view and time 
of operation, respectively, and dVjdz is the differential comov¬ 
ing volume. Here, z(Lmax, ‘5) is the maximum redshift at which a 
burst with would still be above the instrumental flux limit 
S ; Liim(5, z) is the minimum observable luminosity as a function 
of z (i.e. that corresponds to a flux above S). 

If fiL) and f{z) are known, it is possible to derive from Eq.[T] 
the flux distribution of the population of GRBs observable by a 
given detector, knowing its instrumental parameters. By revers¬ 
ing this argument, one can assume the functional forms of f{L) 
and fiz) (e.g. specified through a set of free parameters) and con¬ 
strain them by fitting the the model flux distribution (i.e. A(> S )) 
to the observed flux distribution of a given instrument. This indi¬ 
rect method has been used to infer the luminosity function (e.g. 
Firmani et al. 2004; Salvaterra & Chincarini 2007; Salvaterra et 
al. 2009, 2012) by fitting e.g. the flux distribution of the large 
population of GRBs detected by BATSE. 

The number of free parameters, if both fiL) and f{z) are to 
be constrained, can be large. One possibility is to assume that, 
based on the massive star progenitor origin of long bursts, the 
GRB cosmic rate traces the cosmic star formation rate, i.e. f{z) 
oc fi^iz). In this way, the method allows to derive the free param¬ 
eters of (p{L) by fitting the result of Eq. [T] to large, statistically 
significant, samples of observed GRBs. The assumed f{z) can 
be tested by fitting the observed redshift distribution of a sample 
of bursts with measured z. 

In the simplest scenario, the two functions (p(L) and if/(z) are 
independent. However, more realistic analyses also considered 
the possible evolution of either the luminosity function or the 
GRB formation rate with redshift. For example, in the case of 
luminosity evolution, the burst luminosity depends on z through 
the relation L(z) = To(l + z)'^ (luminosity evolution scenario). 
Alternatively, the GRB formation rate if/(z) varies with redshift 


with a similar analytical dependence fiz) oc (/c*(z)(l H- z)‘' (den¬ 
sity evolution scenario). This means that the progenitor charac¬ 
teristics evolve with z and that the ratio of the GRB formation 
rate to the cosmic star formation rate is not constant. Among 
the drawbacks of this method is that it relies on the assumption 
of a specific functional form of (p(L) (and simple functions, e.g. 
power law, broken power law or power law with a cutoff at low 
luminosities, have been adopted) and it allows to test only for 
evolution of the luminosity or of the rate independently. 

Salvaterra et al. (2012) applied the indirect method to a com¬ 
plete sample of GRBs detected by Swift (§3). They And that 
either a luminosity evolution with k = 2.1 + 0.6 or a density 
evolution with d = 1.7 + 0.5 can reproduce the flux distribu¬ 
tion of BATSE bursts and the redshift distribution of the Swift 
complete sample. However, they can not discriminate between 
these two scenarios. They derive the luminosity function 0(L) 
testing two analytical models: a power-law with an exponen¬ 
tial cut-off at low luminosities and a broken power-law (BPL, 
fiL) cx (LILff’^, where a and b are the slopes of the power 
law below and above the break Lb) adopting a minimum GRB 
luminosity Lmin = 10"'^ erg/s. For the BPL model they found 
a = -0.741};36, b = -1.921°;]^, Lb = 5.5^®;^ x 10^° erg/s and 
a = b = 2.32:';°'^’, Lb = 3.8+“ x 10 ^^ g^g/s in the case 

of luminosity and density evolution scenario, respectively. 

The alternative method is based on the direct derivation of 
the 0(L) and tf/{z) from observed samples of GRBs with mea¬ 
sured z and L. This method has been inherited from the studies 
of the luminosity function of quasars and blazars (e.g. Chiang 
et al. 1998, Maloney & Petrosian 1999, Singal et al. 2012,2013) 
and it has been applied to GRBs (Lloyd et al. 1999, Kocevski & 
bang 2006). Wanderman & Piran (2010) adopt a maximum like¬ 
lihood estimator to derive the discrete luminosity function and 
cosmic formation rate. They use the sample of ~100 GRBs de¬ 
tected by Swift with measured redshift (through optical afterglow 
absorption lines and photometry). Despite this sample might suf¬ 
fer from incompleteness, they derive fiL) (extending from 10^" 
erg/s up to 10^"' erg/s) which can be represented by a broken 
power law with a = -1.2, b - -2.4 and Lb = 10^^ ^ erg/s. Sim¬ 
ilarly, they also derived the discrete GRB formation rate fiz) 
which can be represented by a broken power-law as a function 
of (1 -i- z) with indices ni = 2.1+°^ and n 2 - -L4+j g peaking at 
z = 3.1+gg. This rate is consistent with the SFR of Bouwenset al. 
(2009) forz < 3. They assume that the luminosity is independent 
from redshift. 

More recently, Yu et al. (2015 - Y15 hereafter) and Pet¬ 
rosian et al. (2015 - P15 hereafter) applied a statistical method 
to reconstruct the discrete (f>{L) and fiz) from a sample of Swift 
bursts with measured redshifts. They both And a strong lumi¬ 
nosity evolution with k ~ 2.3. Their results converge towards 
a cumulative luminosity function described by a broken power- 
law with a = -0.14 + 0.02, f = -0.7 + 0.03, Lb = 1.43 x 10^' 
erg/s (Y15) and a = -0.5, f = -2.2, Lb = 10^' erg/s (P15). 
These indices (a and J3) are the slopes of the cumulative lumi¬ 
nosity function which is linked to the differential one through 

the integral <I>(L) = N(> L) - f{L) dL. Therefore, for a 
BPL luminosity function, the slopes of the differential form are 
(a, b) — (a - l,p - 1). 

Intriguingly, they And, at odd with respect to previous works, 
that the GRB rate is flat or decreases from the local Universe up 
to z = 1. If compared to the SFR, this behaviour would imply a 
relative excess of the GRB formation rate with respect to i/^,(z) at 
z <1 (if both are normalised to their respective peaks). They dub 
this behaviour the excess of GRBs at low redshifts. This result 
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is puzzling also because it is completely at odds with the find¬ 
ings of the works based on the properties of GRB host galaxies. 
In fact, Vergani et al. (2014), Perley et al. (2015a, 2015b, 2015c) 
and Kruhler et al. (2015), performed multi-wavelength and spec¬ 
troscopic studies on the properties (stellar masses, luminosities, 
SFR and metallicity) of GRB host galaxies of different complete 
GRB samples and compared them to those of the star-forming 
galaxies selected by galaxy surveys. All their results clearly in¬ 
dicate that at z < 1 only a small fraction of the star-formation 
produces GRBs. 

Both P15 and Y15 apply a statistical method (Efron & Pet¬ 
rosian 1992) to remove the redshift dependence of the luminos¬ 
ity induced by the flux-cut in the selected GRB sample. They 
use GRBs detected by Swift with measured redshifts. However, 
while Yu et al. (2015) work with the bolometric luminosity 
of GRBs, Petrosian et al. (2015) adopt the luminosity in the 
Swift/BAT (15-150 keV) energy band. Y15 use all GRBs de¬ 
tected by Swift with a measured redshift and well constrained 
spectral parameters: despite their relatively large number of ob¬ 
jects (~ 130), this is an incomplete sample. P15 account for in¬ 
completeness by cutting their sample to a relatively large flux 
level, at the expense of the total number of bursts, i.e working 
with ~ 200/250 events with measured z. 

Independently from the method adopted to recover (p(L) and 
i//(z), one key point is the definition of the sample. S12 showed 
the importance of working with complete samples of GRBs (see 
also Hjorth et al. 2012). Here we start with the so called BAT6 
Swift sample (S12) and extend it with additional bursts that sat¬ 
isfy its selection criteria (§3). We will then use it to derive the 
luminosity function and the cosmic GRB formation rate (§4). 


3. BATS extended version 

The BAT6 complete sample as defined in S12 was composed by 
58 Swift GRBs with (i) favourable observing conditions for their 
redshift measurement as proposed in Jackobsson et al. (2006) 
and (ii) a peak photon flux P > 2.6 ph cm“^ s“' (integrated over 
the 15-150 keV Swift/BAJi energy band). This sample, which is 
complete in flux by definition, turned out, after selection, also to 
be highly complete (~90%) in redshift (i.e. 52/58 bursts have z). 

The study of the isotropic equivalent luminosity L of the 
bursts of the BAT6 sample requires the knowledge of their broad 
band prompt emission spectrum. Nava et al. (2012) collected the 
46/52 GRBs, within the BAT6, with measured and z. Six 
bursts with measured z did not have measurements. One of 
the main drawbacks of the narrow/soft energy range of the BAT 
instrument is the difficulty to measure the peak Ep of the vEy 
spectrum for several bursts it detects. Other instruments (e.g. 
Konus/Wind - Aptekar et al. 1995, Eerm//GBM - Meegan et al. 
2009 or SuzakuPf^AM. - Yamaoka et al. 2009) compensate for 
this lack, thanks to their wide energy range, measuring a spec¬ 
trum extending from few keV to several MeV. 

Sakamoto et al. (2009) showed that for Swift bursts with mea¬ 
sured £p there is a correlation between the slope of the spectrum 
apL (when fitted with a single powerlaw model) and the peak en¬ 
ergy £p (measured by fitting a curved model). With the aim of 
enlarging the sample of Nava et al. (2012), we estimated £p of 
six bursts of the BAT6, whose BAT spectrum is fitted by a sin¬ 
gle power law, through this relation (Sakamoto et al. 2009) and 
verified that the values obtained are consistent with those of the 
other bursts (we performed the K-S test finding a probability 
of ~ 70% that the two sets of peak energies originate from the 
same distribution). We find that all but one GRB (i.e. 070306) 


have £p®®* = £“'’*( 1 -H z) consistent with the upper/lower limit re¬ 
ported in Nava et al. (2012). Therefore, we firstly extended the 
BAT6 sample of Nava et al. (2012) with measured z and L to 
50/58 bursts. 

Since the construction of the BAT6, other bursts satisfying its 
selection criteria were detected by Swift. Moreover, some burst 
already present in the original BAT6 sample were re-analysed 
and either their redshifts and/or their spectral properties were 
revised. So our first aim was to revise the BAT6 sample. In par¬ 
ticular, the revision of 8 redshifts is here included (marked in 
italics in the table - their luminosity has been updated). The re¬ 
vised BAT6 sample then contains 56/58 GRBs with measured z 
and 54/58 with also a bolometric isotropic luminosity L. Con¬ 
sidering only the redshift, the sample is ~ 97% complete, while 
if we also require the knowledge of L, the completeness level is 
only slightly smaller (~ 93%). 

Then, we extended the BAT6 revised sample with the new 
events detected since 2012 which satisfy the observability crite¬ 
ria of Jackobsson et al. (2006). The extended sample contains 99 
GRBs. We collected the spectral parameters of the new GRBs in 
the existing literature. For those events (6/41 bursts) with only 
a Swift single power law spectrum, we estimated £p through the 
Sakamoto et al. (2009) relation. The BAT6 extended (BAT6ext, 
hereafter) counts 82/99 GRBs with z ( and 81/99 with z and L. 
Its completeness in redshift is ~ 82%. 

The BAT6ext is presented in Tab. [T] The first 58 bursts are 
the original BAT6, while the others constitute the extension. For 
each GRB, Tab. [T] shows the redshift z, the spectral parame¬ 
ters (high and low photon indices a and f and the rest frame 
peak energy the peak flux with the relative energy band, 

the isotropic equivalent luminosity L. The spectrum is a cut-off 
power-law (CPF) if only the low energy photon index a is re¬ 
ported and a Band function if also the high energy photon index 
yS is given. When z is not measured, we report the observed peak 
energy. The luminosities reported in the table are calculated in 
the [1 - 10"^] keV rest frame energy range only for those GRBs 
having both z and Ep. 

4. Luminosity function and GRB formation rate 

In this section we will apply the C method as originally pro¬ 
posed by Fynden-Bell et al. (1971) and applied to GRBs by e.g. 
Yonetoku et al. (2004,2014), Kocevski & Liang (2006), Wu et 
al. (2014). This method is based on the assumption that the lu¬ 
minosity is independent from the redshift. However, as discussed 
in Petrosian et al. (2015) a strong luminosity evolution could be 
present in the GRB population. Efron & Petrosian (1992) pro¬ 
posed a non-parametric test to estimate the degree of correlation 
of the luminosity with redshift induced by the flux in a flux- 
limited sample. This is also the case of the BAT6ext sample and 
the first step will be to quantify the degree of correlation. Yu et 
al. (2015) and Petrosian et al. (2015) indeed found that the lumi¬ 
nosity evolves with redshift within their samples as (1 -H z)^ ‘*^-»3* 
(Y15) or (1 +2)2-3 ±o.8 (P 15 ), We applied the same method of Y15 
and P15 (also used in Yonetoku et al. 2004,2014) to the BAT6ext 
sample: we define the luminosity evolution L - Lo(l + zi^ (as in 
Y15), where Lo is the de-evolved luminosity, and compute the 
modified Kendall correlation coefficient (as defined in Efron & 
Petrosian 1992). We find, consistently with the results of Y15 
and P15, k - 2.5. Similar results were obtained through the in¬ 
direct method (see §2) by S12 using the BAT6 sample. 

We can now define the de-evolved luminosities Lo = T/(l + 
z)^ for every GRBs and apply the Lynden-Bell C method to 
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Table 1. BAT6ext (BAT6 extended) GRB complete sample. Columns report, in order, the redshift z, the spectral photon indices a and /?, the 
peak flux in units of 10“’ erg cm“’ s“* (except for those with the * which are in units of photons cm“’ s“‘), the respective energy hand, the 
rest-frame peak energy and the holometric equivalent isotropic luminosity L (calculated in the [1 - 10“*] rest-frame energy range). For L we 
also give the Icr error. *’ Bursts with missing z are reported with their observer frame for completeness, they are not used in the present 
work. The peak energy has been estimated with the relation of Sakamoto et al. (2009). The eight GRBs (already present in the compilation of 
Nava et al. 2012) with an updated redshift estimate are marked in italics. In the last column we report the references, in order, for the spectral 
parameters and for the redshift: 1) Nava et al. (2012), 2) Covino et al. (2013), 3) Kruhler et al. (2015), 4) GCN# 12133, 5) GCN# 12135-12137, 6) 
GCN#12190, 7)GCN#12352, 8)GCN#12424, 9)GCN#12431, 10) GCN# 12749, 11) GCN# 12761, 12) GCN#12801, 13) GCN# 12839, 14) 
GCN# 12874, 15) GCN# 12865, 16) GCN# 13120, 17) GCN# 13118, 18) GCN# 13412, 19) GCN# 13536, 20) GCN# 13532, 21) GCN# 13559, 
22) GCN# 13562, 23) GCN# 13634, 24) GCN# 13628, 25) GCN# 13721, 26) GCN# 13723, 27) GCN# 13990, 28) GCN# 13992, 29) 
GCN# 13997, 30) GCN# 14052, 31) GCN# 14419, 32) GCN# 14437, 33) GCN# 14487, 34) GCN# 14491, 35) GCN# 14469, 36) GCN# 14493, 
37) GCN# 14545, 38) GCN# 14575, 39) GCN# 14567, 40) GCN# 14720, 41) GCN# 14808, 42) GCN# 14796, 43) GCN# 14869, 44) 
GCN# 14959, 45) GCN# 14956, 46) GCN# 15064, 47) GCN# 15145, 48) GCN# 15144, 49) GCN# 15203, 50) GCN# 15187, 51) GCN# 15413, 
52) GCN# 15407, 53) GCN# 15452, 54) GCN# 15450, 55) GCN# 15669, 56) GCN# 15805, 57) GCN# 15800, 58) GCN# 15853, 59) 
GCN# 16134, 60) GCN# 16125, 61) GCN# 16220, 62) GCN# 16217, 63) GCN# 16262, 64) GCN# 16310, 65) GCN# 16423, 66) GCN# 16473, 
67) GCN# 16495, 68) GCN# 16489, 69) GCN# 16512, 70) GCN# 16505. 


GRB 

z 

a \ fi ] 

Peak flux 

Range 


L 

Ref 




10 “^ erg cm “^ 

* phot cm “^ 

( keV ) 

( keV ) 

( xlO ^^ erg / s ) 


050318 

1.44 

- 1.34 ± 0.32 

2.20 ± 0.17 

[15 - 150 ] 

115 ±27 

4.76 ± 0.86 

1,1 

050401 

2.9 

- 1 . 0 [- 2 . 45 ] 

24.5 ± 1.2 

[20 - 2000 ] 

499 ± 117 

201 ± 11 

1,1 

050416 A 

0.653 

- 1 . 0 [- 3 . 4 ] 

5.0 ± 0 . 5 * 

[ 15 - 150 ] 

26 ±4 

0.97 ± 0.12 

1,1 

050525 A 

0.606 

- 0.99 ± 0.11 

47.7 ± 1 . 2 * 

[15 - 350 ] 

127 ±6 

7.24 ± 0.28 

1,1 

050802 ' 

1.71 

- 1 . 6 ± 0.1 

2.21 ± 0.35 

[ 15 - 150 ] 

301 

9.51 ± 1.71 

1,1 

050922 C 

2.198 

- 0.83 ± 0.24 

45 ±7 

[20 - 2000 ] 

416 ± 118 

187 ±30 

1,1 

060206 

4.048 

- 1 . 12 ± 0.30 

2.02 ± 0.13 

[ 15 - 150 ] 

409 ± 116 

49.6 ± 7.1 

1,1 

060210 

3.91 

- 1 . 12 ± 0.26 

2.8 ± 0 . 3 * 

[ 15 - 150 ] 

574 ± 187 

52 . S ± Il.l 

1,1 

060306 

1.55 

- 1 . 2 ± 0.5 

4.71 ± 0.28 

[15 - 150 ] 

178.5 ± 76.5 

11.49 ± 2.26 

1,2 

060614 

0.125 

- 1.5 

11 . 6 ± 0 . 7 * 

[ 15 - 150 ] 

55 ±45 

0.05 ± 0.01 

1,1 

060814 

1.92 

- 1.43 ± 0.16 

21.3 ± 3.5 

[ 20 - 1000 ] 

750 ± 245 

71 . 7 ± 13.1 

1,1 

060904 A 

- 

- 1.22 ± 0.05 

13 ±3 

[ 20 - 10000 ] 

235 ± 25 '’ 

- 

1,1 

060908 

1.88 

- 0.93 ± 0.25 

2.81 ± 0.23 

[15 - 150 ] 

426 ± 207 

12 . 7 ± 3.1 

1,1 

060912 A " 

0.94 

- 1.85 ± 0.08 

25 ±9 

[ 20 - 10000 ] 

127 

20.6 ± 7.4 

1,1 

060927 

5.47 

- 0.81 ± 0.36 

2.47 ± 0.17 

[ 15 - 150 ] 

459 ± 90 

108.7 ± 13.1 

1,1 

061007 

1.261 

- 0.75 ± 0.02 [- 2.79 ± 0 . 09 ] 

120 ± 10 

[ 20 - 10000 ] 

965 ± 27 

109.2 ± 8.9 

1,1 

061021 

0.346 

- 1 . 22 ± 0.13 

37.2 ± 9.3 

[20 - 2000 ] 

1046 ±485 

1.77 ± 0.46 

1,1 

061121 

1.314 

- 1.32 ± 0.05 

128 ± 17 

[20 - 5000 ] 

1402±185 

142 ± 19 

1,1 

061222 A 

2.09 

- 1.00 ± 0.05 [- 2.32 ± 0 . 38 ] 

48 ± 13 

[ 20 - 10000 ] 

1091 ± 167 

140 ± 38 

1,1 

070306 " 

1.50 

- 1.67 ± 0.1 

3.04 ± 0.16 

[15 - 150 ] 

> 263 

> 9.99 

1,1 

070328 

2.063 

- 1.11 ± 0.04 [- 2.33 ± 0 . 24 ] 

59 ± 12 

[ 20 - 10000 ] 

2349 

157.6 ± 37.6 

1,3 

070521 

2.087 

- 0.93 ± 0.12 

41 . 2 ± 9.1 

[20 - 1000 ] 

685.6 ± 73.6 

144.1 ± 32.6 

1,3 

071020 

2.145 

- 0.65 ± 0.29 

60.4 ± 20.8 

[20 - 2000 ] 

1013 ±204 

213±73 

1,1 

071112 C " 

0.82 

- 1.09 ± 0.07 

8.0 ± 1 . 0 * 

[ 15 - 150 ] 

596 

6.57 ± 0.86 

1,1 

071117 

1.331 

- 1 . 53 ± 0.15 

66.6 ± 18.3 

[20 - 1000 ] 

648 ±317 

95.4 ± 28.4 

1,1 

080319 B 

0.937 

- 0.86 ± 0.01 [- 3.59 ± 0 . 45 ] 

226 ±21 

[20 - 7000 ] 

1307 ±43 

101 . 6 ± 9.4 

1,1 

080319 C 

1.95 

- 1.20 ± 0.10 

33.5 ± 7.4 

[20 - 4000 ] 

1752 ± 504 

96.1 ± 21.7 

1,1 

080413 B 

1.10 

- 1.23 ± 0.25 

14.0 ± 0.6 

[15 - 150 ] 

163 ± 34 

14 . 9 ± 1.8 

1,1 

080430 " 

0.77 

- 1.73 ± 0.08 

1.82 ± 0.13 

[ 15 - 150 ] 

149 

1 . 16 ± 0.13 

1,1 

080602 ^ 

1.820 

- 0.96 ± 0.63 

19.2 ± 5.8 

[ 20 - 1000 ] 

1216 

51 ± 17 

1,3 

080603 B 

2.69 

- 1.23 ± 0.64 

15.1 ± 3.9 

[ 20 - 1000 ] 

376 ±214 

116.6 ± 38.9 

1,1 

080605 

1.64 

- 1.03 ± 0.07 

160 ±33 

[20 - 2000 ] 

665 ± 48 

308.7 ± 62.8 

1,1 

080607 

3.036 

- 1.08 ± 0.06 

269 ± 54 

[20 - 4000 ] 

1691 ±169 

2260 ± 446 

1,1 

080613 B 

- 

- 1.05 ± 0.18 

47.6 ± 13.1 

[20 - 3000 ] 

33 ± 239 '’ 

- 

1,1 

080721 

2.591 

- 0.96 ± 0.07 [- 2.42 ± 0 . 29 ] 

211 ±35 

[20 - 7000 ] 

1785 ±223 

1039 ± 173 

1,1 

080804 

2.20 

- 0.72 ± 0.04 

7.30 ± 0.88 

[8 - 35 000 ] 

810 ±45 

27.0 ± 3.2 

1,1 

080916 A 

0.689 

- 0.99 ± 0.05 

4.87 ± 0.27 

[8 - 35 000 ] 

208 ± 11 

1.08 ± 0.06 

1,1 

081007 

0.53 

- 1.4 ± 0.4 

2.2 ± 0 . 2 * 

[25 - 900 ] 

61 ± 15 

0.43 ± 0.09 

1,1 

081121 

2.512 

- 0.46 ± 0.08 [- 2.19 ± 0 . 07 ] 

51.7 ± 8.3 

[8 - 35 000 ] 

608 ± 42 

195.4 ± 33.7 

1,1 

081203 A 

2.10 

- 1 . 29 ± 0.14 

2.9 ± 0 . 2 * 

[15 - 150 ] 

1541 ±756 

28.3 ± 8.9 

1,1 

081221 

2.26 

- 0.83 ± 0.01 

24.2 ± 0.5 

[ 8 - 35000 ] 

284 ±2 

101 ± 2 

1,3 

081222 

2.77 

- 0.90 ± 0.03 [- 2.33 ± 0 . 10 ] 

17.6 ± 0.58 

[8 - 35000 ] 

630 ±31 

95 ±6 

1,1 

090102 

1.547 

- 0.97 ± 0.01 

29.3 ± 0.91 

[8 - 35000 ] 

1174 ±38 

45.7 ± 1.4 

1,1 

090201 

2.1 

- 0.97 ± 0.09 [- 2.80 ± 0 . 52 ] 

73.0 ± 12.6 

[20 - 2000 ] 

489.8 

269.5 ± 47.6 

1,2 

090424 

0.544 

- 1.02 ± 0.01 [- 3.26 ± 0 . 18 ] 

91 . 2 ± 1.4 

[8 - 35000 ] 

250.0 ± 3.4 

11 . 16 ± 0.18 

1,1 

090709A 

1.8 

- 0.85 ± 0.08 [- 2.7 ± 0 . 24 ] 

39 ±6 

[20 - 3000 ] 

834.4 

91 . 9 ± 13.9 

1,2 

090715 B 

3.00 

- 1.1 ± 0.37 

9.0 ± 2.5 

[20 - 2000 ] 

536 ± 164 

82.6 ± 25.2 

1,1 

090812 

2.452 

- 1.03 ± 0.07 

2.77 ± 0 . 28 * 

[ 100 - 1000 ] 

2023 ± 663 

96.3 ± 16.0 

1,1 

090926 B 

1.24 

- 0.19 ± 0.06 

4.73 ± 0.28 

[ 8 - 35000 ] 

212.0 ± 4.3 

4.28 ± 0.25 

1,1 

091018 

0.971 

- 1.53 ± 0.48 

4.32 ± 0.95 

[ 20 - 1000 ] 

55 ±26 

4.75 ± 1.33 

1,1 

091020 

1.71 

- 1.20 ± 0.06 [- 2.29 ± 0 . 18 ] 

18 . 8 ± 2.6 

[ 8 - 35000 ] 

507 ± 68 

32.7 ± 5.2 

1,1 

091127 

0.49 

- 1.25 ± 0.05 [- 2.22 ± 0 . 01 ] 

93.8 ± 2.3 

[8 - 35000 ] 

51 . 0 ± 1.5 

9.09 ± 0.24 

1,1 

091208 B 

1.063 

- 1.29 ± 0.04 

25.6 ± 0.97 

[8 - 35000 ] 

246 ± 15 

17.5 ± 0.7 

1,1 

100615A 

1.4 

- 1.24 ± 0.07 [- 2.27 ± 0 . 11 ] 

8.3 ± 0 . 2 * 

[8 - 1000 ] 

206.4 ± 20.4 

10.15 ± 0.87 

1,2 

100621 A 

0.542 

- 1.70 ± 0.13 [- 2.45 ± 0 . 15 ] 

17 . 0 ± 1.3 

[20 - 2000 ] 

146 ± 23 

3.17 ± 0.34 

1,1 

100728 B 

2.106 

- 0.90 ± 0.07 

5.43 ± 0.35 

[ 8 - 35000 ] 

404 ± 29 

18.7 ± 1.3 

1,1 

110205 A 

2.22 

- 1.52 ± 0.14 

5.1 ± 0.7 

[20 - 1200 ] 

715 ±238 

25.1 ± 4.3 

1,1 

110503 A 

1.613 

- 0.98 ± 0.08 [- 2.7 ± 0 . 3 ] 

100 ± 10 

[20 - 5000 ] 

572 ±50 

180.7 ± 19.7 

1,1 

110709 A 

- 

- 1 . 16 ± 0.02 

15 . 4 ± 1 . 7 * 

[ 10 - 1000 ] 

533 ± 37 '’ 

- 

4 ,- 

110709 B 

<4 


11 ± 1 

[20 - 5000 ] 

278 ):"'’ 

- 

5,6 

110915 A 

- 

- 0.94 ± 0.23 

3.3 ± 0 . 2 * 

[15 - 150 ] 

124.8 ± 41 . 4 '’ 

- 

7 ,- 
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GRB 

z 

a[/S] 

Peak flux 

10“^ erg cm“^ 

* phot cm“^ 

Range 

(keV) 

£'p 

(keV) 

L 

(xlO^* erg/s) 

Ref 

111008A" 

4.989 

-1.86 ±0.09 

6.4 ± 0.7* 

[15 - 150] 

384 

303.3 ±48.8 

8,9 

111228A" 

0.715 

-1 [-2.27 ±0.06] 

12.4 ±0.5* 

[15-150] 

46 

3.64 ± 0.27 

10,11 

120102A 

- 

-1.19 ±0.03 

22.8 ± 1.6* 

[10- 1000] 

380 ± 33" 


12- 

120116A 

- 

-1.31 ±0.41 

4.1 ±0.3* 

[15-150] 

- 

- 

13,- 

120119A 

1.728 

-0.98 ±0.03 [-2.36 ±0.09] 

16.9 ±0.4* 

[10- 1000] 

189.2 ±8.3 

56.9 ±2.7 

14,15 

120326A 

1.798 

-1.41 ±0.34 

4.6 ± 0.2* 

[15-150] 

115± 19 

10.8± 1.8 

16,17 

120703A 

- 

-o.si!“ 

54± 13 

[20- 10000] 

295+88b 

- 

18,- 

120729A" 

0.8 

-1.62 ±0.08 

2.9 ±0.2* 

[15-150] 

192 

1.27±0.11 

19,20 

120802A 

3.796 

-1.21 ±0.47 

3.0 ± 0.2* 

[15-150] 

274.3 ± 93.0 

40.7 ±5.7 

21,22 

120811C 

2.671 

-1.4±0.3 

4.1 ±0.2* 

[15-150] 

157.5 ±20.9 

25.4 ±4.5 

23,24 

120907A 

0.97 

-0.75 ± 0.25 

4.3 ± 0.4* 

[10- 1000] 

304.4 ±64.8 

2.45 ± 0.27 

25,26 

121123A 

2.7 

-0.96 ±0.2 

2.6 ± 0.2 

[15-150] 

240.5 ± 18.9 

14.9 ± 1.4 

27,28 

121125A 

- 

-1.38 ±0.06 

4.2 ±0.3 

[10- 1000] 

196 ± 26" 

- 

29,- 

121209A" 

2.1 

-1.43 ±0.08 

3.4 ±0.3* 

[15-150] 

494 

17.4± 1.6 

30,3 

130420A 

1.297 

-1.52 ±0.25 

3.4 ±0.2* 

[15-150] 

76.3 ±15.6 

3.77 ± 0.65 

31,32 

130427A 

0.339 

-0.958 ±0.006 [-4.17 ±0.16] 

6900 ± 100 

[20- 1200] 

1371.3 ± 10.7 

384.1 ±5.7 

33,34 

130427B‘' 

2.78 

-1.64±0.15 

3.0 ±0.4 

[15-150] 

386 

29 ±5 

35,36 

130502A 

- 

-1.0 ±0.3 

7± 1 

[8- 1000] 

83 ± 17" 

- 

37,- 

130505A 

2.27 

-0.69 ±0.04 [-2.03 ±0.03] 

690 ± 30 

[20 - 1200] 

2063.4 ± 101.4 

3959 ± 172 

38,39 

130527A 

- 

1.04 ±0.04 

500 ± 30 

[20- 10000] 

1380±120" 

- 

40,- 

130606A" 

5.913 

-1.14±0.15 

2.6 ±0.2* 

[15-150] 

2032 

229.5 ± 24.3 

41,42 

130609B 

- 

-0.66 ±0.22 [-2.6 ±0.2] 

13.6 ±0.4* 

[10- 1000] 

491 ±20" 

- 

43,- 

130701A 

1.555 

-0.9 ±0.21 

17.1 ±0.7 

[15-150] 

227.9 ±31.6 

28.9 ± 1.5 

44,45 

130803A 

- 

0.85 ±0.09 

7.1 ±0.3* 

[10- 1000] 

141.6 ± 12.2" 

- 

46,- 

130831A 

0.479 

-1.61 ±0.06 [-3.3 ±0.3] 

25 ±3 

[20- 10000] 

81.3±5.9 

3.68 ± 0.45 

47,48 

130907A 

1.238 

-0.91 ± 0.02 [-2.34 ± 0.07] 

220 ± 10 

[20- 10000] 

881.8 ±24.6 

185.6 ±8.8 

49,50 

131030A 

1.293 

-0.71 ±0.12 [-2.95 ±0.28] 

100 ± 10 

[20- 10000] 

405.9 ±22.9 

103 ± 11 

51,52 

131105A 

1.686 

-0.88 ±0.38 [-2.33 ±0.33] 

20 ±2 

[20- 10000] 

419 ± 102 

37.4 ±4.9 

53,54 

140102A 

- 

-0.71 ±0.02 [-2.49 ±0.07] 

49.7 ±0.5* 

[10- 1000] 

186 ±5" 

- 

55,- 

140206A 

2.73 

-1.04 ±0.15 

19.4 ±0.5* 

[15-150] 

376.4 ±54.1 

141.5 ±4.8 

56,57 

140215A 

- 

-0.66 ±0.11 [-2.94 ±0.35] 

35.7 ±3.5 

[20- 10000] 

214 ± 14" 

- 

58,- 

140419A 

3.956 

-0.63+Sg [-2.3+“] 

47+18 

[20- 10000] 

1452.1 ±416.3 

572.6 ±25.2 

59,60 

140506A 

0.889 

-0.9 ±0.2 [-2.0 ±0.1] 

14.2 ±0.7* 

[10- 1000] 

266 ± 68 

11.5± 1.3 

61,62 

140512A 

0.72 

-1.33 ±0.03 

11.0±0.3* 

[10- 1000] 

1011±145 

5.28 ± 0.47 

63,64 

140619A 

- 

-1.45 ±0.14 

4.6 ±0.2* 

[15-150] 

117.8 ±46.2" 

- 

65,- 

140628A 

- 

-1.56 ±0.09 

2.8 ± 0.2* 

[15-150] 

- 

- 

66,- 

140629A 

2.275 

-1.42 ±0.54 

4.7 ± 0.7 

[20- 10000] 

281.7 ±57.4 

27.1 ±5.5 

67,68 

140703A 

3.14 

-1.10±0.06 

4.1 ±0.2* 

[10- 1000] 

732.8 ±58.0 

41.6±2.2 

69,70 




l+z Lohrg/s] 


Fig. 1. Left panel: GRB formation rate obtained with the C method using the BAThext sample (black solid line). The dashed red line and 
the dot-dashed cyan line are the SFR models of Hopkins & Beacom (2006) and Cole et al. (2001) shown here for reference. All the curves are 
normalised to their maxima. Right panel, luminosity function 0(Lo) obtained with the method using the BAT6ext sample (black solid line). 
The best fit model describing this function is a broken power-law (dashed green line) with {a = -1.32 ± 0.21, b = -1.84 ± 0.24, Li, = 
erg/s). The orange dot-dashed line is the luminosity function obtained by S12 in the case of pure luminosity evolution. 


derive the cumulative luminosity function 0(Lo) and the GRB 
formation rate i/^(z). 

For the ith GRB in the BAT6ext sample, de¬ 
scribed by its (Lo,i,2i), we consider the subsample 
J\ - {y'lf'O.j > f-o,! A Zj < Zmax,i} and call N\ the number 
of GRBs it contains. Similarly, we define the subsample 
J[ = [j\Lo.i > Liim.i A Zj < Zi} and we call Mi the number of 


GRBs it contains. Lum i is the minimum luminosity correspond¬ 
ing to the flux limit S of the sample at the redshift Zi- Zmax.i is the 
maximum redshift at which the i-th GRBs with luminosity ^ 
can be observed (i.e. with flux above the limit S). Lum and Zmax 
are computed applying the K-correction because the limiting 
flux is computed in the observer frame 5w//f/BAT [15 - 150] 
keV energy band. 
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Through Mi and Ni we can compute the cumulative luminos¬ 
ity function <l>(Lo) and the cumulative GRB redshift distribution 

dz)- 




J<1 


( 2 ) 


and 



(3) 


From the latter we can derive the GRB formation rate as: 


l/'iz) = —t—(1 +2) 
dz 

where dV(z)ldz is the differential comoving volume. The 
differential luminosity function 0(Lo) is obtained deriving the 
cumulative one <l>(Lo)- 

The functions 0(Lo) and (//{z) are shown in Fig. [1] Errors on 
0(Lo) are computed propagating the errors on the cumulative one 
assuming Poisson statistics. The errors on the ij/iz) are computed 
from the number n of GRBs within the redshift bin. We assume 
that the relative error e - II ^Jn is the same affecting 


dV(z) 


dz 


(4) 


5. Results 

The luminosity function <p{Lo) obtained with the BAT6ext sam¬ 
ple is shown in Fig. [T] (right panel) by the black symbols. 
Data are normalised to the maximum. The best fit model (green 
dashed line in Fig. [1] - right panel) is represented by a broken 
power law function with a - -1.32 + 0.21, b = -1.84 + 0.24, 
Lb = jq51.45±o.15 erg/s (where a and b represent the slopes of 
the power law below and above Lb - ;if^/d.o.f. = 0.47). In Fig. 
[T]we also show, for comparison, the luminosity function derived 
in S12 through the indirect method described in §2 in the case 
of pure luminosity evolution (a - b - 

Lb = 5.5;|;3 4 X 10^° erg/s). This model is consistent with the re¬ 
sult obtained in our analysis. 

The GRB formation rate i/r(z) obtained with the BAThext is 
shown by the black symbols in Fig. [T] (left panel). The green 
dashed and the cyan dot-dashed lines are, respectively, the SFR 
of Hopkins & Beacom (2006) and Cole (2001). Data and models 
are normalised to their peak. Contrary to the results reported in 
P15 and Y15, the ^{z) that we derive increases up to z ~ 2. 
This trend is consistent with the independent estimates obtained 
through host galaxy studies. 

The BAT6ext sample has a smaller completeness in redshift 
(~ 82%) with respect to the revised BAT6 (~ 97%). We checked 
if this could in some way modify the shape of ^(z). For this 
reason we also computed the \}/{z) using only the 56 objects of 
the revised BAT6 sample which turns out to be slightly steeper 
both at low and high redshifts than the one obtained with the 
BAT6ext but totally consistent within the errors. We conclude 
that the lower completeness in redshift of the BAT6ext does not 
introduce any strong bias in the i/f(z) and <^(Lo) obtained. 


6. Monte Carlo test of the C~ method 

We now test the method used to derive i/'(z) and (p{Lo). 
Through a Monte Carlo simulation (similar to e.g. Ghirlanda et 
al. 2015) we explore how well the method adopted in ^can 
recover the input assumptions, i.e. ^{z) and LF <p{L). In partic¬ 
ular, we will show that if the sample used is highly incomplete 
the resulting GRBFR and LF can differ significantly from the in¬ 
put ones. In particular, incomplete samples (either in flux and/or 
redshift) may produce a misleading excess of low redshift GRBs 
with respect to the assumed ^{z). 

We simulate GRBs distributed in redshift according to the 
GRB formation rate ipiz) of Li 2008 (see also Hopkins & Bea¬ 
com 2008): 


0.0157-H 0.118z 
1 -H (z/3.23)4-66 


(5) 


where if/iz), in units of Mq yr“' Mpc“^, represents the formation 
rate of GRBs and we assume it can extend to z < 10. We stress 
that for the scope of the present test any other functional form of 
i/'(z) could be assumed. 

We adopt a luminosity function (piL) as obtained by Sal- 
vaterra et al. 2012 from a complete sample of Swift GRBs: 


(p{L) oc 



, L < Lb 
, L > Lb 


(6) 


composed by two power-laws with a break at Lb. We adopted ar¬ 
bitrary parameter values: a = -1.2, L = -1.92 and Lb = 5 x 10^° 
erg s '. We further assume an evolution of the luminosity propor¬ 
tional to (1 H- z)^ with k = 2.2. We assume this is the bolometric 
luminosity of the simulated bursts and compute the correspond¬ 
ing bolometric flux. Also for cpiL) we use this functional form 
but any other function could be assumed for the scope of the 
present test. 

With these two assumptions we simulate a sample of GRBs 
with a flux limit Liim = 5 x 10“^ erg cm“^ s“* and we analyse it 
with the C method. Accounting for the truncation of the sample 
we recover through the statistical method of Efron & Petrosian 
(1992) the luminosity evolution in the form (1 -Hz)*, with k ~ 2.2. 

Then we work with the de-evolved GRB luminosities Lq = 
L/(l + z)^'^ and derive the GRB formation rate ip{z) and the 
luminosity function (piLo) through the C method proposed by 
Lynden-Bell et al. (1971). The left panel of Lig.l^shows (black 
symbols) that we recover the GRB formation rate of Eq. |5] 
(dashed green line) that we adopted in the simulation. Similarly 
the right panel of Eig. |3 shows that we also recover the lumi¬ 
nosity function that we adopted in the simulated sample (Eq.|^- 
shown by the dashed green line in Eig.|2]i. 

We then tested what happens if we apply the same method 
to an incomplete sample. Eirstly we applied the C method to 
the same simulated sample (which is built to be complete to 
Liim = 5 X 10“* erg cm“^ s“*) from which we removed ran¬ 
domly a fraction of the bursts close to Liim. This new sample is 
clearly incomplete to Fn^. The results are shown in Eig.|2]by the 
blue symbols. We find that the GRB formation rate ip{z) is flat 
at low redshifts (i.e. below z = 2) showing a clear excess with 
respect to the assumed function (cf blue symbols with the green 
dashed line in the left panel of Eig.|2]l. The luminosity function 
is flatter than the assumed one (cf blue symbols with the green 
dashed line in Eig.|2]right panel). 

Similar results are obtained by assuming for the derivation 
of i/r(z) and 0(Lo) a flux limit which is a factor of 5 smaller than 
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Fig. 2. Left: GRB formation rate (normalised to its peak) for the simulated population of GRBs with flux limit 5 x 10“* erg cm“^ s“* (black 
symbols). The GRB formation rate assumed in the simulation is shown by the dashed green line. The red symbols show the results obtained from 
the same sample using, for the analysis, a flux limit a factor of 5 smaller than the real one. Blue symbols are obtained by mimicking the sample 
incompleteness by removing randomly some GRBs near the flux threshold adopted for the sample selection. Right panel: cumulative luminosity 
function, normalised to the first bin. The black, red and blue symbols are the same as for the left panel. The assumed luminosity function is shown 
by the dashed green line. 


that used to construct the simulated sample. This is another way 
to make the sample artificially incomplete. The results are shown 
by the red symbols in the panels of Fig. |2] Note that in this sec¬ 
ond test the sample used is the same but it is analysed through 
the C method assuming it is complete with respect to a flux 
limit which is smaller (a factor of 5) than the one corresponding 
to its real completeness (i.e. 5 x 10“* erg cm“^ s“*). 

These simulations show that if the samples adopted are 
highly incomplete in flux, one obtains an excess at the low red- 
shift end of the GRB formation rate and a flatter luminosity func¬ 
tion. 

7. Summary and discussion 

The aim of this work is to derive the luminosity function of long 
GRBs and their formation rate. To this aim we apply a direct 
method (Lynden-Bell et al. 1971) and its specific version already 
applied to GRBs, e.g. Yonetoku et al. (2004, 2014), Kocevski & 
Liang (2006), Wu et al. (2014), P15, Y15. This is the first time 
this method is applied to a well defined sample of GRBs which 
is complete in flux and 82% complete in redshift. 

We build our sample of long GRBs starting from the BAT6 
complete sample (Salvaterra et al. 2012): this was composed by 
58 GRBs detected by the Swift satellite satisfying the multiple 
observational selection criteria of Jackobsson (2006) and having 
a peak photon flux P > 2.6 ph cm“^ s *. Here, we update the 
redshift measurement of 8 GRBs of the BAT6 (marked in ital¬ 
ics in Tab. [T]i and accordingly revise their luminosities. Then, 
we update this sample to GRB 140703A ending with 99 objects. 
We collect their redshift measurements and spectral parameters 
from the literature (see Tab.[T]i. The BAT6ext sample has a red¬ 
shift completeness of ~ 82% (82/99 burst with z measured) and 
counts 81/99 bursts with well determined L. 

We analyzed the BAT6ext sample searching for a possible 
luminosity evolution induced by the flux threshold through the 
method proposed by Efron & Petrosian (1992). We find that the 
L - z corTelation introduced by the truncation due to the flux 
limit can be described as L = L()(l + z)^ with k = 2.5. This result 
is in agreement with what found by other authors (Yonetoku et 


al. 2004, Wu et al. 2012, P15, Y15). With the BAT6ext sample 
after de-evolving the luminosities for their redshift dependence 
we find that: 

- the luminosity function fiLo) is a monotonic decreasing 
function well described by a broken powerlaw with slopes 
a = -1.32 + 0.21 and/? = -1.84 +0.24below and above,re¬ 
spectively, a characteristic break luminosity Lb = 

erg/s (green dashed line in the right panel of Fig.[TJ. This re¬ 
sult (shape, slopes and characteristic break) is consistent with 
the luminosity function found by S12 (orange dot-dashed 
line in Fig.[T]- right panel). 

- The cosmological GRB formation rate if(z) (black solid line 
in the left panel of Fig. [T]) increases from low redshifts to 
higher values peaking at z ~ 2 and decreases at higher red- 
shifts. This trend is consistent with the shape of the SFR 
of Hopkins & Beacom (2006) and Cole et al. (2001) (color 
lines in Fig. [T]|. Our results on ij/(z) is in contrast with the 
GRBFR recently found by P15 and Y15, who, applying the 
same method to differently selected GRB samples, report the 
existence of an excess of low redshift GRBs. 

However, the result that GRBs evolve in luminosity is not the 
probe that GRBs had experienced a pure luminosity evolution. In 
fact the C method assumes the independence between L and z 
and the non parametric method of Efron & Petrosian used to get 
the de-evolved luminosities assigns all the evolution in the lu¬ 
minosity. For this reason we are not able to distinguish between 
an evolution in luminosity or in density (see also Salvaterra et 
al. 2012). The possible density evolution case requires the in¬ 
vestigation of the applicability of similar methods to the GRB 
samples and will be the subject of a forthcoming work (Pescalli 
et al. , in preparation). 

Finally, intrigued by the different results with respect to Y15 
and P15, we performed Monte Carlo simulations in order to test 
the robustness of the C method. We showed that this method 
can correctly recover the LF 4>{Lq) and the GRBFR f{z) assumed 
in the simulation only if the sample of GRBs it is applied to is 
complete in flux and has a high level of completeness in red¬ 
shift. Using incomplete samples or a sample incomplete in red- 
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shift, the resulting GRBFR and LF can be different from the as- [] 
sumed ones. Indeed, this could account for the excess of the rate 11 
of GRBs at low redshift as recently reported. 
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Fig. A.l. The blue filled points represent the observed CLF S(L) while 
the orange squares represent the true CLF t(L). The solid blue line (a = 
-1.39 ± 0.14, b = -2.38 ± 0.41, U = lO^^'^*^ erg/s) and the dashed 
orange line (a = -1.22 ± 0.1, b = -2.09 ± 0.95, Lb = 10^^'^*' erg/s) are 
the best fit models of the observed and true CLF respectively. 


Appendix A: The redshift integrated iuminosity 
function 

In this section we show how we derive the luminosity function of 
GRBs integrated over all the redshift space using the BAThext 
sample. This is not the typical luminosity function which is de¬ 
rived through indirect methods, it is free from any functional 
form and only uses the l/Vmax concept, i.e. the maximum co¬ 
moving volume within which a real burst with an observed lumi¬ 
nosity could be detected by a given instrument. This is a gener¬ 
alisation of the < y/Vniax > method proposed by Schmidt (1968) 
and applied to quasars (see Avni & Bachall 1980 for an exhaus¬ 
tive description). 

However, the luminosity function obtained with this method 
is not the canonical (f>{L) but it is the result of the integration of 
the latter convolved with the GRB formation rate (//{z) over z. We 
call this function the convolved luminosity function S(L) (CLF). 
In principle the luminosity function could evolve with redshift 
{(p{L, z)). For this reason the shape of the CLF could be different 
from the shape of the canonical (p{L). It would be simply pro¬ 
portional to (p(L) only if the luminosity function does not evolve, 
either in luminosity or in density. Indeed, it is possible to express 
2(L) in terms of luminosity function and GRB formation rate; 


2(L) = 



0(L, z)(A(z) dz 


(A.l) 


The advantage is that it can be obtained directly from the 
data and it is extremely robust if derived through a flux limited 
sample. Any evolution with redshift of the luminosity or density 
is in this CFL. We use the 81/99 GRB with both z and determined 
L reported in Tab. 1. 

For each GRB of the BAThext sample with an associated L 
we estimate the maximum volume Vmax within which the burst 
could still be detected because its flux would be larger than the 
chosen threshold i.e. Pn^ = 2.6 ph cm“^ s ' in the 15-150 keV 
energy band. The observed photon flux in the Swift/BAOC energy 
band as a function of the varying redshift is; 


P(z) = 




And^iz) 


L 


\l(\+z)keV 


EN{E) dE 


(A.2) 


where N{E) is the observed photon spectrum of each GRB and 
difz) is the luminosity distance at redshift z. The extremes of 
the integral in the denominator correspond to the same values 
adopted to compute L. The maximum redshift Zmax corresponds 
to the redshift satisfying F’(zmax) = Fum. 

Considering the typical Swift held of view G = 1.33 
sterad and the time of activity of Swift T ~ 9 yr that covers the 
BAT6ext sample, we can compute the rate pi = AnlQ.TViaax,i for 
each GRB. We divided the observed range of luminosities in bins 
of equal logarithmic width A and estimate £j(Lj) for each bin as 

= (A.3) 


where the sum is made over the bursts with luminosities Lj - 
A/2 < Li < Lj-i-A/2. The discrete convolved luminosity function 
is shown in Fig. lA.ll The normalisation is obtained also consid¬ 
ering that the bursts in the BAT6ext sample represent approxi¬ 
mately 1/3 of the total number of Swift detected GRBs with peak 
flux P>2.6 ph cm“^ s"'. We verified that the BAT6ext sample (99 
objects) is representative in terms of peak flux distribution of the 
larger population. The error bar associated to the discrete CLF 
are mainly related to the Poissonian error on the objects count¬ 
ing within the luminosity bin (see also e.g. Wolter et al. 1994). 

When computing individual rates, we have used the mission 
elapsed time T. However, this is the observer frame time and the 
rate should be corrected for the cosmological time dilation. This 
means that, at higher z, the same subset of sources should occur 
with a larger frequency. Therefore, we average out the elapsed 
time on the redshift interval [0, Zmax] 


< r >= 


Fornax T 7 

Jo 

Jo 


(A.4) 


and use it in the computation of Eq. IA.3I As shown in Fig. lA.ll 
the true CLF £(L) is flatter than the observed one because, in 
general, the true elapsed time is smaller than the observed one. 
Moreover, the correction is more pronounced for high luminosity 
GRBs observable up to high redshifts. 

The CLF obtained with the extended BAT6ext sample is 
shown in Fig. lA.ll bv the filled symbols. The observed CFL (blue 
symbols in Fig. lA.ll i can be adequately represented by a bro¬ 
ken power law function (solid blue line in Fig. lA.ll) with slopes 
-1.39 + 0.14 and -2.38 + 0.41 below and above the break lumi¬ 
nosity Lb = erg/s (^^/d.o.f. = 1.04). When we correct 

for the cosmological time dilation the true CLF appears slightly 
flatter (slopes -1.22 + 0.1 and -2.09 + 0.95 below and above the 
break Lb = 10^^'^="' erg/s - ;tVd.o.f. = 0.99). 

We can look at this function as a redshift integrated luminos¬ 
ity distribution. This is the most direct information that we can 
obtain from data, in fact, it is obtained without any type of as¬ 
sumption or observational constraints. Similarly to the flux dis¬ 
tribution log A - log 5 and to the redshift distribution it can be 
used as a constrain. The LF and the GRBFR obtained with other 
methods should reproduce this CLF once convolved together and 
integrated over z. 
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